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We present preliminary results for the chiral behavior of charged pseudo-Goldstone-boson masses 
and decay constants. These are obtained in simulations with Nf=2-\-\ flavors of tree-level, 0{a)- 
improved Wilson sea quarks. In these simulations, mesons are composed of either valence quarks 
discretized in the same way as the sea quarks (unitary simulations) or of overlap valence quarks 
(mixed-action simulations). We find that the chiral behavior of the pseudoscalar meson masses in 
the mixed-action calculations cannot be explained with continuum, partially-quenched chiral per- 
turbation theory. We show that the inclusion of 0{a^) unitarity violations in the chiral expansion 
resolves this discrepancy and that the size of the unitarity violations required are consistent with 
those which we observe in the zero-momentum, scalar-isotriplet-meson propagator. 
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1. Introduction 

The objective of our collaboration is to calculate hadronic observables which are relevant for 
determining fundamental quark properties, such as quark masses or quark flavor-mixing and CP- 
violation parameters, and to do so with controlled extrapolations to the physical limit of Nf =2+1 
flavor QCD, where Mji ~ 135 MeV, the lattice spacing a vanishes and the volume is infinite. To 
achieve that goal we consider two approaches. In both, the seas are composed of Nf =2+1 flavors 
of tree-level, 0(a) -improved Wilson (W) quarks [|l]]. We perform "unitary" simulations where the 
valence quarks are discretized in the same way as the sea quarks. We also perform "mixed-action" 
calculations, with overlap Ginsparg -Wilson (GW) [^] valence quarks. In the latter, the 

valence sector possesses a full, continuum-like chiral symmetry which greatly simplifies the 
renormalization of electroweak operators, such as those encountered in neutral kaon mixing. It 
also guarantees that matrix elements are automatically 0(<3)-improved, to the extent that the sea 
quarks are. As discussed by Stefan Krieg at this conference [^, recent advances have allowed us 
to perform Nf =2+1 simulations, for instance, down to Mji ~ 190MeV with a ~ 0.09 fm and in 
cubic volumes of side L ~ 4.2 fm. Thus, we expect to be able to reach the near-continuum chiral 
p-regime of Gasser and Leutwyler without the conceptual problems of staggered fermions [^, ^. 
This means that we should be able to extrapolate lattice results to the physical point in a model- 
independent way, by using Wilson [|l^-|l^, partially quenched (PQ) chiral perturbation theory 
(ZPT) [0-0] for the unitary simulations, and mixed-action (MA) PQ;tPT [|l8[ [Bl, [l^ for the 
GW-on-W simulations. 

One of the drawbacks of using a mixed-action approach is the presence of discretization- 
induced unitarity violations. Fortunately, it should be possible to account for the low-energy man- 
ifestations of these violations with MAPQ;t;PT. We present here preliminary results for the quark- 
mass dependence of the masses and leptonic decay constants of the pseudo-Goldstone bosons 
(PGBs) of chiral symmetry breaking. In particular, we investigate the effects of unitarity violations 
in these quantities, as obtained in our mixed-action simulations, and attempt to correlate these ef- 
fects with those which we observe in the scalar, isotriplet, propagator, where they are expected 
to be particularly large [pQ|]. Of course, our study of the masses and decay constants of the PGBs 
is primarily motivated by the very interesting phenomenology they give rise to. They allow the 
determination of a variety of fundamental quantities, such as light quark masses, the ratio of CKM 



elements |V„i/Kj| [21] and important LO and higher-order low-energy constants (LECs) of the 
effective chiral Lagrangian. However, we postpone the presentation of results for these quantities 
to later pubUcations. 



2. Finite-volume mixed-action PQ;tPT and unitarity violations 

As shown in | |20| ] for the case of Nf degenerate flavors, the propagator of the ao is affected by 
potentially large 0{a^) unitarity violations in a mixed-action scenario. For the case which interests 
us, with Nf =2+1, to simplify expressions we suppose that the light sea {£) and valence (v) quark 
masses are tuned such that masses of the corresponding PGBs are equal: M^, = Ma = Mj^. We 
denote the strange sea quark by s, and by M„ the mass of charge pseudoscalar mesons composed 
of two valence quarks with the mass of the s. ' Then, LO MAPQ;^PT gives for the zero-momentum, 

'From now on we call these mesons "non-singlet ss" pseudoscalar mesons. 
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ao propagator: 

Qo(?,A) = Zj{gQ,aX)a^Y,i^2qiix,t)qiq2iO)) (2.1) 

X 

t^+o. B^(X) r , , 2 , , a^A , 

where the valence quarks are chosen to be degenerate, with m\ = m2 = m^, and where we have 
assumed, for simplicity, that the time and space extent of the lattice are infinite. The quantity B is 
the condensate over the square of the decay constant F. Zs{go,aX) is a renormalization constant 
for scalar densities and A is a QCD renormalization scale. The functions Cxy(0 denote the zero- 
momentum propagators of the two-particle states, XY , and are given by Cxr(?) = exp[— (M^ + 
My)f] / {4MxMy) at LO. There are two physical contributions, coming from intermediate two-kaon 
and 717] states. The unitarity violations are of order and are proportional to a quantity A, which 
has mass dimension four. These violations only vanish in the continuum limit. Moreover, they are 
exponentially and polynomially enhanced in t, also compared to the contribution from the oq not 
shown in Eq. At asymptotic times, the unitarity violations are the dominant contribution. 

Unitarity violations also affect PGB masses and decay constants, but only at NLO. Let us 
consider a pseudoscalar meson composed of two distinct quarks with masses mi and m2. Then, 
according to NLO MAPQ;^PT, the square of this meson mass has the following generic form: 

(M22)gLO = (mi+m2)B{l + (4^[pQ-logs(At,Mn,M22,M«,M,,) 

+ (2a6 - 0:4) (m ) (2M2^ +Ml) + (2a8 - as ) (m )Mj2 + ^ (2.2) 

+aj8M + a2Ax{UV-logs(At,Mii,M22) + 7M(M)} } , 



where B and A are defined after Eq. (2.1); "PQ-logs" and "UV-logs" denote partially-quenched 



and unitarity-violating quenched-like logarithms, respectively; /X is the renormalization scale in 



the effective theory and, in Eq. (2.2), the quark masses and the B must both be either renormal- 
ized in the same scheme in QCD or bare; 7m(m) is an a priori unknown counter-term; the LECs 
a,- are related to the original Gasser and Leutwyler constants thi^ough a;(/i) = 8(47r)^L,(jU); FV 
stands for finite- volume corrections; Pm is a mass-dimension three quantity which parameterizes 
0{a) discretization errors and whose parametric size will be specified below. Note that in fits to 
lattice results obtained at a single lattice spacing, the discretization errors proportional to aj^M and 
a^A7M(jU) get absorbed into the LO LEG, B. 

In applying the general form of Eq. (^^), it is useful to distinguish three cases of interest. 
In the continuum or in the case of GW valence on GW sea quarks, mi and m2 can be taken to 
be the Lagrangian masses and A and jSm identically vanish, i.e. there are no unitarity violating 
nor 0{a) discretization errors. In our W-on-W unitary simulations, we take mi and m2 to be 
the "measured" axial Ward identity (AWI) masses. Moreover, the constant j^M is 0(Aq(.q) if the 
fermions used are straight Wilson fermions, C?(asAQ^p) if they are tree-level 0(a) -improved as 
they are in our simulations, or zero if they are non-perturbatively 0(<3)-improved. Finally, A = 
and discretization-induced unitarity violations are absent. In the mixed-action, GW-on-W case, mi 
and m2 can be taken to be the valence, overlap Lagrangian masses. 15m has the same parametric 
size as for the unitary simulations, depending on the level of improvement of the sea. However, in 
the mixed-action case, A does not vanish a priori. Thus, we expect that meson masses will suffer 
from discretization-induced unitarity violations at finite a. 
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The NLO expression for the decay constant, F12, of a charged pseudoscalar meson is simi- 
lar in form to Eq. (2.2), with the factor {nii +m2)B replaced by F and where the "PQ-logs"and 
corresponding counter-terms and finite-volume corrections are changed to the partially-quenched 
expressions appropriate for decay constants. The corresponding general form applies to the three 
cases discussed above in much the same way, with the substitutions jSm jSf and Jm ^ 7f- The 
main difference appears in the GW-on-W case, where the mixed-action unitarity violations propor- 
tional to a^A are purely valence SU (3) -flavor breaking, and do not depend on jJ.. 



3. Results for the uq propagator and the charged PGB masses and decay constants 

The results discussed below were obtained from simulations which were described at this con- 
ference by Stefan Kiieg [^. We recapitulate here only the main ingredients, omitting algorithmic 
considerations. The gauge action used is tree-level Symanzik improved [^]. The sea quarks are 



described by six-step stout-smeared []23|], tree-level C?(a)-improved Wilson fermions For the 
unitary, W-on-W simulations, the valence quarks are discretized in the same way. In the GW-on-W, 
mixed-action case, valence quarks are three-step HYP-smeared [ p4| ] overlap fermions [|2|- J, with 
a negative mass parameter p = 1 . 

We have performed five, 2+1 sea-flavor simulations at a lattice spacing a ~ 0.09 fm (j8 = 
3.57). In these simulations, the mass of the charged pions composed of two light sea quarks are 
Mjc ~ 190, 300, 410, 490 and 570MeV. To keep finite-volume errors small, all simulations are 
performed in cubic three-volumes with sides L such that M^^L > 4. The strange quark mass used 
in our simulations is slightly overestimated: with such a strange quark, the mass of a kaon, when 
extrapolated in light-quark mass to the physical point, is approximately 7% higher than the physical 
kaon mass. There are 34 gauge configurations at ~ 190 MeV, 68 at Mji ~ 300 MeV and 0(100) 
at the three other simulation points. In the mixed-action case, the overlap quark masses are chosen 
such that the mesons which they compose are approximately degenerate with those composed of 
the corresponding Wilson sea quarks. At Mj^ ~ 190 and 300 MeV, we have a second overlap valence 
strange quark whose mass is approximately 30% smaller than that of the strange sea quark. 

We begin with preliminary results for the zero-momentum, ao propagator, Cag(f), defined in 



Eq. (2J_). In Fig. |I] we plot the unrenormalized propagators, C^^^{t), as a function of Euclidean 



time t, obtained in the two GW-on-W simulations with the lightest u and d quarks. The propagators 
go negative at relatively short times and then asymptotically go back up to zero. Moreover, the 
effect is less pronounced for the simulation with the more massive quarks, a trend which persists 
as one increases the u and d quark masses further. This behavior is qualitatively consistent with 



the prediction of MAPQ;^PT given in Eq. (2.1), assuming A > 0. The agreement can be made 



quantitative also. To verify this, we perform fits of the propagators, at asymptotic times t, to the PQ 
and FV generalization of Eq. ( p] ) for the bare Cao{t). The fits have only one parameter, namely 
a^A. For the pre-factor {B/ZsY, we take Mj2/(?wi +m2), where M12 is the "measured" meson 
mass and mi ^2 are the bare masses of the GW quarks which compose it. The kaon, pion and rj 
masses which appear in the expression for C^^^{t) are constrained to take on the values obtained 
from prior fits to kaon and pion two-point functions. These mesons are composed of a sea and 
a valence quark. We obtain their masses by combining the corresponding valence-valence and 
sea-sea masses at LO in the chiral expansion. This means that the parameter a^A that we fit here 
contains contributions from an 0{a^), mixed-action operator in the chiral Lagrangian. Fortunately, 
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Figure 1: Bare, zero-momentum propagator of the ao as a function of time over half the time extent of 
our lattices, as obtained in our two GW-on-W simulations with the lightest u and d quarks. The solid curves 



represent our best fits to the partially-quenched and finite-volume generalization of Eq. ( 2. 1 ) in the fit region, 
and the dashed curves their extensions to earlier times. 



if the same prescription is used for obtaining valence-sea meson masses which appear at NLO in 
quantities such as or ^n, the unitarity violations there will be parameterized by the same a^A. 

The one-parameter fits to C^^^{t) for Mj^ ~ 190 and 300 MeV are performed for t/a in the 
range [12,32], where t/a = 32 is the midpoint of our lattices in both cases. The results of these fits 
are plotted in Fig. |l|. As the figure suggests, both fits have good x^/dof. The values obtained for 
the unitarity-violation parameter are a^A = 0.015(6) and 0.024(10), respectively, for M^j ~ 190 and 
300 GeV, and are thus consistent. For the lattice spacing at which the simulations are performed, 
these values correspond to aVA ~ 0.27 GeV and 0.35 GeV. Since aVA competes with pion, kaon 
and T] masses in the chiral expressions for and F12 in the mixed-action case, it is clear that 
these unitarity- violating contributions cannot be neglected a priori. 

We now turn to an analysis of the GW-on-W decay constant. We begin with this quantity, 
because we use the extrapolated aTv to determine the lattice spacing as well as to normalize cor- 
rections in chiral expressions with factors of {4naFji)^. aFu is obtained from the pseudoscalar 
two-point function using the AWI. Thanks to the chiral symmetry of the overlap, valence quarks, 
no renormalization is required. This is a simple example of the simplifications brought about by 
the use of a mixed action with chirally symmetric valence quarks. 

There are 21 lattice points for aFi2, of which 5 correspond to charged "pions", 7 to "kaons" 
and 9 to "non-singlet 55" pseudoscalar mesons. We fit these results to the NLO chiral expression 
described in Sec. ^ So as to remain, as much as possible, within in the range of applicability of 
NLO ;^PT, we include in the fit only the 4 lightest pion and 4 lightest kaon points, with Mji < 
500MeV and Mk < 590MeV. The fit has four parameters, which are F, a4(M^), as^Mj^) and a^A. 
Since the ao propagator is more sensitive to unitarity violations than are the decay constants, we fix 
a^A to the value 0.024(10), through a Gaussian prior in the x^- The NLO expression describes the 
data well. Moreover, the resulting value of a^A is 0.025(8), confirming that the chiral behavior of 
the decay constants is consistent with the presence of unitarity violations of the size observed in the 
ao propagator. The value of aFj^ obtained from a self-consistent extrapolation to the physical point, 
using the infinite-volume, continuum fitted function, yields a lattice spacing of a = 0.088(1) fm. 
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Figure 2: The bare condensate ratio as a function of squared sea-pion mass, in lattice units, for the unitary 
simulations (left) and for the mixed-action simulations (right). The vertical scales in the two plots are 
only equal up to a ratio of renormalization constants. The three sets of data (circles) in the W-on-W case 
correspond, from top to bottom, to pion-like, kaon-like and non-singlet .s.f-like pseudoscalar mesons. In the 
GW-on-W plot, there are one additional kaon and two extra ss points at the two smallest M^f's, corresponding 



to the additional valence strange quark that we consider in these simulations. The fits to Eq. (2.2) are plotted 
as line-segments around each fitted point. The physical curves are obtained from the fits by removing the 
FV effects, and in the mixed-action case, the partial-quenching effects and unitarity-violating logs. The pion 
curves are obtained by setting Mii=M22=M(( and M^s—Mf^^^ , the latter being the physical, non-singlet ss 
pseudoscalar meson mass; the kaon curves, by setting Mi i =M(£ and M22=Mss=Mf^^^; the SU{3) curves, by 
setting Mii~M22=Mss—M(('. The vertical dotted lines mark, from left to right, the chiral limit, the physical 
pion and the physical kaon points. 



where the error is statistical. 

Next we consider chiral fits to the unitary W-on-W results for the condensate ratio aB^f^ = 
{aMxjf' / {am\ +<3m2)^^j. The lattice results for this quantity are plotted in Fig. ^ as a function 
of the squared sea-pion mass, M|^. The different sets of points correspond to "pions", "kaons" 
and "non-singlet ss" pseudoscalar mesons. We fit these results to the NLO expression of Eq. (2.2) 
with A = 0. We include in the fit only the six points with < 500MeV and Mk < 590 MeV. 
Here, the fit has only three parameters, namely B, {la^ — 054) (M^) and {la^ — a5)(M,^). The NLO 
expression describes the chiral behavior of aB\f^ very well. 

We now turn to the condensate ratio obtained in the mixed-action, GW-on-W simulations. 
Here, aB^f^ = {aM\2Y / {cim\ +am2f'^^^, where M\2 is the valence meson mass and m\^2 the 
corresponding bare overlap Lagrangian masses. The results for this quantity are plotted in Fig. ^ 
again as a function of M|^. As the plot indicates, the behavior of aB\f'^ here deviates significantly 
from that obtained in the unitary case. Moreover, some of the features of this behavior, such as 
the large increase of aB\2^ for the "pion" points at small Mf^, cannot be explained with only 
continuum PQ chiral logarithms: a divergent term at small appears to be required. Fortunately, 
such a contribution is provided by the unitarity violations discussed in Sec. ^ and exhibited in 
Eq. (2.2). We thus fit the lattice results to the NLO expression of Eq. (2.2), including the unitarity 
violating term proportional to a^\. As for the decay constant fit, we fix fl^A to the value 0.024(10), 
obtained from the propagator, through a Gaussian prior in the x^. There are four parameters 
in the fit, one more than in the W-on-W case. These are B, {la^ — a4){Mn), {2ag, — a5)(Mf,) 
and a constrained a"* A. Again, only the eight points with Mji < 500 MeV and Mk < 590 MeV are 
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included. The description of the condensate ratio given by our NLO chiral expression is good. The 
value of a'^A returned by the fit is 0.020(6). As already noted after Eq. (2.2), in a fit performed 



at fixed lattice spacing, the discretization error proportional to 7m (m) gets absorbed into B, which 
thereby acquires a spurious ;^PT /i -dependence. Of course, this dependence will be eliminated, 
along with all other discretization errors, when the renormalized values of B, obtained at different 
lattice spacings, are extrapolated to the continuum limit. The value of /i chosen here is Mr). A 
lower value will raise the physical curves whereas a larger one will lower them. It is worth noting 
that this spurious -dependence cancels at NLO in ratios such as ms/niud or (w)%=2/('?'?)a'/=3- 
Moreover, it only affects the fitted LECs (2a6 — a^){Mrf) and {2a^ — a5)(M^) very mildly. We 
find very good agreement between the GW-on-W and W-on-W results for these quantities. 

These observations, together with the other results reported on here, suggest that unitarity 
violations are present in our mixed-action results, and that we can subtract them with MAPQ;^PT. 
Of course, simulations at other lattice spacings are required to confirm this conclusion. 
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